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The simulation of advanced solidification processes via digital computer techniques has gained 
widespread acceptance during the last decade or so. Models today can predict transient temperature 
fields, fluid flow fields, important microstructural parameters, and potential defects in castings. 
However, the lack of accurate thermophysical property data on important industrial alloys threatens to 
limit the ability of manufacturers to fully capitalize on the technology's benefits. This paper describes 
a study of the sensitivity of one such numerical model of a steel plate casting to imposed variations in 
the data utilized for the thermal conductivity, specific heat, density, and heat of fusion. The sensitivity 
of the data's variability is characterized by its effects on the net solidification time of various points 
along the centerline of the plate casting. Recommendations for property measurements are given and 
the implications of data uncertainty for modelers are discussed. 

Simulation results can only be as good as the boundary conditions/heat transfer coefficients 
applied and the thermophysical properties used, i.e., density, specific heat, latent heat, and thermal 
conductivity. Several other authors have investigated the important effects of interfacial heat transfer 
coefficients and found a strong dependency of the evolution of the thermal field on these coefficients. 
However, few modelers cite more than a passing concern for the lack of thermophysical data and the 
uncertainty inherent in the data that is available. Depending on the source of thermophysical property 
data, considerable errors may be present. For example, the data for the thermal conductivity of 
tungsten reported in the literature exhibits as much as 300% scatter even though many data sets are 
reported to be accurate to within 1% or less! In addition, the data for a great many complex alloys of 
industrial interest are simply unavailable. 

The casting chosen for this study is the solidification of a horizontal steel plate in green sand 
with a riser at one end as shown in Figure 1. This is the same geometry investigated by Minakawa et 
al l in their study of the feeding of porosity. The nominal properties used were identical to those of 
Minakawa et al and are shown in Table I. The thickness of the plate is 25 mm and the the length is 200 
mm. Centerline shrinkage has been observed to occur in the region indicated in Figure l. 2 The 
assumptions made by Minakawa et al were repeated here. (1) The mold is instanteously filled with 
molten metal at the pouring temperature. (2) The thermal contact resistance at the metal-sand mold 
interface is negligible. (3) Segregation is neglected. (4) The latent heat is released uniformly between 
the liquidus temperature and the solidus temperature irrespective of cooling rate effects. 

The model was tested by comparing calculations of the plate centerline temperatures with the 
measured temperatures of Bishop and Pellini. 2 Excellent agreement was found in the temperature 
range between the solidus and the liquidus. 

The effects of various uncertainties upon the predicted solidification times due to uncertainties 
in the steel's thermal conductivity, specific heat, latent heat and density are shown in Figures 2, 3, 4, 
and 5 respectively. A comparison of these effects is presented in Figure 6 where the solidification 
times of the node at 57.1 mm from the plate edge are shown. Increasing the thermal conductivity value 
of the steel by 100% only decreased the time to solidify by about 13%. An error of +50% in the 
specific heat of the steel would cause an increase in the predicted time to freeze of -25%, whereas a 
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+50% error in the latent heat would increase the predicted solidification time by 55%. The sensitivity 
of the model is even greater to variations in the density used. The predicted solidification time 
increased by 90% with a +50% error in the steel's density. 

In terms of the molten steel's properties, the model investigated was most sensitive to 
uncertainties in the steel's density and least sensitive to it’s thermal conductivity. The approximate 
sensitivity coefficients are: 
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Although high temperature molten alloys are experimentally difficult, many techniques exist for 
determining most of these properties required to accuracies of the order of ±5% or better. However, 
there is little agreement on stoadanl techniques, little publicly available data on most common industrial 
alloys, and apparently little incentive for researchers to worry about the absolute validity of the data 
they use. Inaccurate data leads to inaccurate results and can stop a development program in its tracks. 
All investigators should strive to critically assess the amount of uncertainty in the data that they use and 
to quantify the expected effects of that uncertainty in their results. The maturation of computer 
modeling from a research tool to a design tool demands no less. 

Conclusions 

1 . The sensitivity of computer solidification models to uncertainties in thermophysical properties can 
be readily assessed by straightforward variation of their values in a one-at-a-time manner. 

2 . Large errors in some of the input thermophysical properties can lead to corresponding large errors 
in computer models' predictions. 

3 . For variations up to ±25% in the thermophysical properties, the relative importance of the input 
properties is: density stee ] > latent heat^i > specific heat^^i > thermal conductivity Steel- 
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Figure 2 

Effects of Uncertainty in Thermal Conductivity 
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Calculated solidification time, t f , min Calculated solidification time, t f , min 


Figure 3 

Effects of Uncertainty in Specific Heat 



Figure 4 

Effects of Uncertainty in Latent Heat 
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Figure 5 

Effects of Uncertainty in Density 









Heat conduction is governed by the following well known 
unsteady equation: 
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pc p 


where T is the temperature 
p is the density 
cp is specific heat 
Lf is the latent heat 
t is time, 

K is the thermal conductivity, and 
the properties are independent of temp. 


A simple-minded, one-at-a-time variation of the properties 
in a validated model would quickly give valuable insight 
into which thermophysical properties are the dominant 
ones governing heat flow in casting. 


40 





Vol. 3, Supp. No. 1 , American Chenical Society and American 
Institute of Physics, 1974, pp. 1-369. 


41 



.150 



Table I 

Nominal Values of the 
Physical Data Used in the Calculations 


Property or Parameter Units Steel Mold 


Specific heat, Cp 
Density of liquid, p 
Density of solid, p 
Thermal conductivity, k 
Heat transfer coefficient, H s 
Latent heat, Lf 
Liquidus temperature, T L 
Solidus temperature, T s 
Pouring temperature, Tp 
Emissivity, e 


J/kg«K 

840 

1050 

kg/m 3 

7100 


kg/*n 3 

7500 

1650 

W/m°K 

31 

1.55 

W/m 2 °K 


20.9 

JAg 

2.7 X 10 5 


°K 

1780 


°K 

1736 


°K 

1868 


- 

0.45 



S. Mlnakawa, I.V. Samarasekera, and F. Weinberg, MetaM. Trana. B 
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Plate Castin 


TEMPERATURE 
CONTOUR PLOT 


LEGEND 

— 0 . 1736E+04 

— 0 . 1758E+04 



MINIMUM 

0 . 29300E+03 
MAXIMUM 

0 . 17938E+04 



TIME 0 . 600E+03 


XMIN 0 . OOOE+OO 
XMAX 0 . 630E+00 
YMIN 0. 000E+00 
YMAX 0 . 300E+00 


FIDAP 6.02 
30 Jul 92 
13:02:28 


Plate Castin 



TEMPERATURE 
CONTOUR PLOT 


LEGEND 


— 0 . 1736E+04 

— 0 . 1758E+04 



MINIMUM 

0 . 29300E+03 
MAXIMUM 

0 . 17804E+04 


TIME 0 . 900E+03 



SCREEN LIMITS 

XMIN 0. 000E+00 
XMAX 0 . 630E+00 
YMIN 0. 000E+00 
YMAX 0 . 300E+00 

Y 

L 

FIDAP 6.02 
30 Jul 92 
12:58:14 
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Effects of Uncertainty in Specific Heat 
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Effects of Uncertainty In Thermal Conductivity 
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Calculated solidification time, t ^ min Calculated solidification time, t f , min 


Effects of Uncertainty in Latent Heat 



Distance from plate end, cm 


Effects of Uncertainty in Density 
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Effects of Errors in Molten Steel Properties 



Percentage error In steel property 


The finite element model calculates an effective specific 
heat from a supplied enthalpy vs. temperature curve. 


dH/dt 
“ dT/dt 


Thus Cp is determined at each integration point or nodal 
point: 


H(T„) - H(T„.i) 
Cp ~ T n - T n .i 
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Change in Enthalpy with Different C p 's 
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Effects of Errors in Sand Properties 
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Conclusions 


1 . The effects of uncertainties in input data can be 
quickly assessed by a straightforward (but 
computer intensive) approach. 

2. The calculated cooling rate of the casting model 
is affected by uncertainties in the thermophysical 
data in the following order of influence: 

1 . density 

2. latent neat 

3. thermal conductivity (sand) 

4. specific heat 

5. density and specific heat (sand) 

6. thermal conductivity 

3. Similar analyses will be useful to benchmark the 
importance of properties in fluid flow analyses. 
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